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Electronic Processes at the Breakdown of the Quantum Hall Effect 
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Microscopic processes giving the energy gain and loss of a two-dimensional electron system 
in long-range potential fluctuations are studied theoretically at the breakdown of the quantum 
Hall effect in the case of even-integer filling factors. The Coulomb scattering within a broadened 
Landau level is proposed to give the gain, while the phonon scattering to give the loss. The 
energy balance equation shows that the electron temperature T e and the diagonal conductivity 
a xx exhibit a bistability above the lower critical electric field E c i. Calculated values of E c i as 
well as T e and a xx at E c i are in agreement with the observed values in their orders of magnitude. 
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In the quantum Hall effectBB (QHE), the diagonal 
conductivity a xx is vanishingly small in the low-current 
regime. When the current is increased up to a criti- 
cal value, a xx increases by several orders of magnitude 
within^riarrow range of the current and the QHE breaks 
downE'Crcr Ll spite of extensive studies on the break- 
down of QHE.tr the mechanism of the breakdown has not 
been fully understood. In many samples, a xx increases 
discontinuously at the critical current and exhibits a hys- 
teresis as a function of the current, which suggests that 
the breakdown of QHE belongs to the nonequilibrium 
phase transition. The phases below and above the tran- 
sition can be considered to be homogeneous in a large 
class of samples, since £he critical current is proportional 
to the sample width.Q'El' In this paper, we study theoret- 
ically the nonequilibrium phase transition between ho- 
mogeneous states at the breakdown of the quantum Hall 
effect in the case of even-intege 

Among a variety of theoriesBB [3 
proposed for the mecha^ism-of the breakdown of QHE, 
the hot-electmn | ^heoryBllj'll3 ) showed the existence of 
the hysteresisErB and reproduced the-observed abrupt 
change_of a xx with the current densityO However, this 
theoryliiP employed an observed result to obtain the en- 
ergy gain j-E = (T X xE 2 (j: current density, E: electric 
field) and jiid not study the microscopic process of a xx . 
StudiesEirliH' were made for the process of a xx by con- 
sidering an inter-Landau-level phonon scattering. The 
calculated cr xx , however, had no hysteresis as a function 
of the electric field and the calculated critical electric 
field was at least one order of magnitude larger than the 
observed one. The microscopic process of the energy gain 
and the energy dissipation remains unsolved. 

In this paper we propose a Coulomb scattering within 
a Landau level in a slowly-varying potential as the dom- 
inant electronic process giving the energy gain. The 
potential in the plane of the two-dimensional electron 
system (2DES) is fluctuating due to ionized donors in 
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the layer several hundred A above 2DES, and its im- 
portance in the bmakdpwn of QHE has been discussed 
in the literature.BE3oli3' Fluctuations are in a scale of 
l v b = 0.1/zm (/ v h is the distance between a potential hill 
and a neighboring valley) and have a standard deviation 
of 20meV before the screening .Elf The screening does not 
completely wash out the fluctuations in strong magnetic 
fields due to the discrete nature of the energy spectrum. 
The screened potential b&s-a. reduced width equal to the 
Landau level separation^' Er and, even in a filling fac- 
tor of 2N(N = 1,2, •■•), electrons (holes) populate in 
the N + 1th (Ath) Landau level (we assume the spin 
degeneracy). Therefore the energy dissipation occurs 
within a Landau level and it is much larger than that 
due to inter-Landau-level scatterings. In such a slowly- 
varying potential in strong magnetic fields, closed orbits 
are formed around its hills and valleys, and their typi- 
.is of the order of / v h and much larger than the 
-j.i.L i en g t h l since Z vh ~ 10/ in B = 5T. A hill 
orbit at the center of the Landau level is in close prox- 
imity to neighboring valley orbits and hoppings between 
these orbits along the electric field are the most domi- 
nant process of the energy gain. We propose a Coulomb 
scattering as the dominant mechanism of a hopping, in 
which an electron hops from a hill (valley) orbit to a 
valley (hill) orbit and, at the same time, one of other 
electrons is excited (relaxed). The dominant process for 
the energy loss, on the other hand, is given by acoustic- 
phonon deformation-potential scatterings, again, within 
the Landau level. In our calculations, (1) we consider 
only the activation transport at the Landau-level center 
and neglect the tunneling transport at the Landau-level 
edge, which limits the quantitative validity of our cal- 
culations to the higher T e branch with the larger a xx . 
We also assume (2) the vanishing lattice temperature: 
Tl = 0, and (3) even-integer filling factors. Our calcu- 
lations show that the electron temperature T e and a xx 
exhibit a bistability above the lower critical electric field 
E c i, giving a hysteresis as a function of E. Calculated 
orders of magnitude of E c \ as well as T e and a xx at E c \ 
agree with the experimental ones. 
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In a steady state, the energy balance equation 
P G {T e ,E)=P L {T ei T L ), 



(1) 



holds for the energy gain Pq and the loss Pl, and Pq is 
related with the diagonal conductivity a xx by 

P G {T e ,E)=j-E = a xx (T e )E 2 . (2) 

We assume, following the previous works, 100) that 
a xx depends on E only through T e (E), which is consis- 
tent with the experiment .E3 The electron temperature 
is obtained as a function of E from the energy balance 
equation. 

The electron distribution function f(e) is given by 

exp[(e - ^)/k B T e \ + 1 

Since the filling factor is an even integer, 2N, and the 
spin splitting is much smaller than ksT e , the chemical 
potential fi is at the middle point between the iVth and 
the iV+lth Landau levels, and the electron distribution 
in the iV+lth Landau level and the hole distribution in 
the iVth Landau level for both spins are described by 
the same function. In the following the zero of energy 
is taken at the center of the N + 1th Landau level and 
/it = —hoj c /2 with oj c the cyclotron frequency. The elec- 
tron (hole) occupation of current-carrying states around 
the N + 2th (N — 1th) Landau- level center, f(huj c ), is 
neglected, which means that we restrict our calculations 
to a lower-Te range of the higher-T e branch. 



Fig. 1. A Coulomb scattering among one hill orbit £3 and three 
valley orbits ei, £2, £4 in a slowly- varying potential V. Afi is the 
chemical potential difference between the hill and the valley. 



The proposed electronic process for the energy gain is 
a Coulomb scattering, as illustrated in Fig. 1, from a hill 
orbit £3 to a valley orbit £4 with an excitation from e\ to 
£2 (we use the energy to label an orbit since we consider 
processes within one hill and one valley). We consider 
the excitation at the closest valley and hill from the hop- 
ping electron. Excitations at larger distances are neg- 
ligible since the Coulomb matrix element, which is ap- 
proximately the interaction between a point charge and a 
dipole, decreases with the distance d as d~ 2 and the tran- 



sition rate as d~ A . The guiding-center distance AX cor- 
responding to the energy difference Ae = £2— £1 = £3—64 
is AX = A£/eEi with Ei the local electric field strength 
due to the slowly-varying potential. The transition rate 
is appreciable only when AX ~ I or smaller, since the 
wavefunction overlap integral S(AX) (along the perpen- 
dicular direction) decreases with AX as 



S(AX) = exp(-AX 2 /4l 2 ) = exp(-A£ 2 /4w 2 ^ 



(4) 



where A£ = A£jhw c and v = leEi/hio c . Inter-Landau- 
level scatterings are neglected here since the correspond- 
ing S(AX) is much smaller: S(AX) = 1CT 11 at v = 0.1 
for a quasi-elastic phonon scattering with eEiAX = huj c 
and S(AX) = 10~ 3 for the most dominant Coulomb 
scattering with Ae = Tiuj c /2 (the transition rate is pro- 
portional to S(AX) 2 for phonon scatterings, and to 
S(AX) 4 for Coulomb scatterings). The change of elec- 
tronic states due to the applied electric field is neglected 
since E c \/Ei ~ 0.1. 



The transition rate 
Coulomb scattering is 



12.34 



of the hill-to-valley 



W 



2tt 



J 34 = -^-|( £ l £ 3 |#c|E2£4)| <5(£i+£ 3 -£2-£ 4 ) (5) 



where H c is the Coulomb interaction: H c = e 2 /e\r a — rt| 
with e the dielectric constant. The matrix element of H c 
is given by 

„2 



(eie 3 |tf c |£ 2 £4> = —S(AXYf L I OTb . 
etvh 



(6) 



Jl is the fraction in length of the orbit £4 in which the 
distance to the orbit £3 is around AX. 7 or b represents the 
dependence on orbital configurations. When the orbits 
£1, £2, and £4 are circles with the common center, 7 or b 
is a function of the radius r\ of orbit £1 as well as AX 
(r4 = Zvh/2), and the square root of the average of I 2 rh 
over r\ and AA" is estimated to be 0.4. 

The energy gain Pq per unit area per unit time is 
related by 



Pa = 



N 2 N C 
2l 2 



:p< 



Gi/i 



(7) 



to the energy gain Pqi in the ith pair of a hill and a 
valley within the TV + 1th Landau level with either spin. 
The coordination number N c is the average number of 
valleys to which an electron hops from a hill. Ni = 4 is 
the number of possible states for each orbit with different 
spins and Landau indices. The energy gain due to a 
single hill-to- valley hopping is £q — eE-Ar with Ar the 
difference between their position vectors. The chemical 
potential is assumed to be constant within each hill and 
each valley and the difference between the hill and the 
valley to be £o- Pa is given by 

p Gi = ■U l -h)h{l-U)W C l 2^£G 

£l ,£2, £31^4 

- £ (l-/l)/2(l-/3)/ 4 ^ 1 C 2! 3 4 £ G (8) 



£l ,£2,£3>£4 



with /„ = /(£„) (n = 1,2,4) and f 3 
We make approximations that f{£\) ^~ 



= / (£3 - so) 
/(e 2 ), f(e 3 - 
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£g) ~ /(E3) - £g/'(e3), and f'(e 3 ) ~ /'(0), by assum- 
ing ksTe ^ Ae ~ vhujc and ksTe £g, which hold 
approximately in the higher T e branch around E c ± if we 
use the value of T e (E c i) and E c \ estimated below. Then 
Pa becomes 



P Gi = £ /l(l-/l)[-/'(0)]^l C 2,34£ 



(9) 



In averaging Pq^ , we change the summation to the inte- 
gration over the energy s. We assume that the density 
of states p(e) is slowly varying so that p{e\) ~ ^(£2) and 
p(e 3 ) ~ p(e 4 ) ~ p(0). With use of (e 2 G ) t = (e£Z vh ) 2 /2, 
we obtain 



with 



P G = A G P 2 /(0) [1-/(0)]// 



P = eEl/huj c , and 



d £ p( £ ) 2 /'(£)- 



(10) 



(11) 



(12) 



p(£) is normalized to be unity when integrated within a 
Landau level and p(e) is a dimensionless density of states: 
p(e) = p(e)huj c . 

As the most dominant process of the energy loss at 
vanishing lattice temperature, we consider the acoustic- 
phonon emission due to the deformation potential. The 
rate Wf 2 of the phonon emission with a transition from 
£1 to £ 2 is 

W i* = X E l^| 2 K £ i |ex P (^-r)| £ 2 )| 2 <5(£ 1 - £ 2 - hu q ) 

q 

(13) 

where q = (q x ,q y ,q z ) and w 9 = ciq are the phonon 
wavevector and angular frequency, respectively, with q 
the group velocity of the longitudinal acoustic mode. 
C q = qD(h/2p m Vcu q ) 1 / 2 with D the deformation po- 
tential, p m the density, and V the volume of the sample. 
The matrix element is calculated to give 

, j 1/ r \ 1 via ( 27 3/2 S(AX) 2 

d^dq y \(s 1 \eM^r)\e 2 )\ = L l (1 + g2/b 2 )3 - 

(14) 

with L p the perimeter of the orbit. We have used the 
Fang-Howard wavefunctionEIP along z (perpendicular to 
the plane): ( (z) = (b 3 /2) 1 / 2 z exp(-bz/2). We here put 
q z = q since g/(q 2 + <7^) 1/2 ~ luj c vjci ~ 3 at P=5T. 
We also assume that orbits are circles with radius 7*0 and 
average the transition rate over rrj. 

The energy loss P L per unit area per unit time is 

^{Pu)i (15) 



Pr. 
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vh 



with 



Pr. 



- ^/l(l-/ 2 )WT 2 -(£l-£2) 

£l>£ 2 



(16) 



Using the approximations used already in the calculation 



of P G : /(£i) ~ f{ £ 2) and p(£i) ~ ^(£2), we obtain 

P L = A L f e I f (17) 
with T e — kBT e /huj c and 



-4, 



4%/5 lvh -(Df3) 2 v 3 I p 



(18) 



TT^ 2 h P 

with /3 = (u /1)(lu c /ujo), ul = H/p m cil 2 , and w = cj/i. 
/ p is defined by 

r- 2 

exp(-x 2 / 2 ) (19) 



Ip ~Jo dX {l + a 2 x 2 f 



with a = vuj c /cib. 

Phonon scatterings (absorptions and emissions) also 
give an electron hopping between a hill and a valley and 
contribute to the energy gain. However, the energy gain 
due to phonon scatterings is shown to be 10 -3 of that 
due to Coulomb scatterings even at T L = T e . 



Fig. 2. (a) The energy gain Pq at three values of the electric 
field E and the loss Pl (both divided by A^If) as a function 
of the electron temperature T e (T e = kBT e /Hui c )- Points of in- 
tersection indicates steady states, (b) T e and cr X x {<?xx /cro! f = 
4/(0) [1 — /(0)] with <To the constant) in steady states as a func- 
tion of E (both have the lower branch of T e = and a xx = 0). 
P G (E) = P L (E) oc T e (E) if the weak T e dependence of I } is 
neglected. 



The energy gain Pq and the loss P L are plotted as a 
function of the electron temperature T e in Fig. 2(a). The 
points of intersection (Pg = Pl) give T e in steady states. 
The number of points of intersection increases from one 
to three as increasing the electric field E through a crit- 
ical value E c i, and a bistability appears above E c \ (one 
of three in the middle corresponds to an unstable state). 
The nonlinear T e dependence of Pg giving this bistabil- 
ity is due to the presence of the activation energy for 
the current-carrying states at the Landau-level center. 
T e = is a stable solution at any E and the upper criti- 



1 



Hiroshi Akera 



cal field is P C 2 — oo, since P G = Pl = at T e = in the 
present approximations of Tl = and no tunneling con- 
duction. T e in steady states as a function of E is given 
by the equation Pq = Pl which is simplified to 

AG f, 2 



-T-E* = g(T e ) 
Ah 



T e [exp(l/2T e ) + exp(-l/2T e ) + 2], 

(20) 

and is plotted with a xx (E) = P G /E 2 in Fig. 2(b). At the 
critical point, f ec = f e (E cl ) = 0.32, and A G E 2 cl /A L 
2:2. 

The estimate of the critical electron temperature is 
free from ambiguities in values of A G and Al and is 
given by T ec = 0.32?Lw c //cb. Unfortunately, there exists 
no direct experimental estimate of T e . Indirect estimates 
have been made from the energy balance equation using 
the observed temperature dependence of a xx at E <C 
E cl and a calculated energy lossME^ Theoretical and 
experimental estimates of T e at E c \ are 

T ec (theory) = 32K [B = 5T), 

T ec (ref JH) = 8K (B = 3.8T), 

T ec (refff) = 15K (P = 7.6T). 

The discrepancy between the theory and the experiments 
is reduced if we consider in the theory the nonzero lattice 
temperature and the smaller activation energy P ac . In 
the present theory P ac =&lj c /2, while E ac = 0.7Uw c /2 
from the observed o xx (T) S3 

The estimate of the lower critical field E c \ depends 
on A L /A G . We use B = 5T, l vh /l = 10, b = 0.03A~\ 



v = 0.1, D=l0eV, p m = 5.3xl0 3 kg/ni , c 4 = 5xl0 3 m/s 
0.1 and < I% Th 



N c = 2, f L 
p(e) =const. Then we obtain 



> l/2 = 



3 „ 

0.4, and assume that 



E cl (theory) = 20V/cm (P = 5T), 

P c i(ref.|) = 65V/cm {B = 4.7T), 

P cl (ref|jJ|) = 40V/cm (P = 3.8T), 

P c i(ref|f) = lOOV/cm (B = 7.6T). 

The order of magnitude agrees between the theory and 
the experiments. Since there are large ambiguities in 
values of l v h/l, N c , and Jl, the discripancy between the 
theory and the experiments is within the limitation of 
accuracy (E c i oc N c and N c takes a smaller value in a 
sparse network of conduction channels between hills_and 
valleys, which is probable in disordered potentialsEy ). 
Note that P c i in the real system is larger than the present 
theoretical estimate, since we have neglected fluctuations 
around the homogeneous steady state. 

The estimate of the energy dissipation at E c \ is given 
by P G (T ec ,E cl ) = P L (T ec ), and depends on A L and //. 
Theoretical and experimental estimates of P G at E c \ are 



P Gc (theory) = 5 x lO^Js^cm- 2 (B = 5T), 

P Gc (ref.|) = 1.5 x lO^Js^cm- 2 (B = 4.7T), 
P Gc (ref.[l||) = 0.8 x lO^Js^cm- 2 (B = 3.8T), 
P Gc (ref||) = 0.6 x lO^s^cm" 2 (P = 7.6T). 
The agreement in P Gc is poorer than in E&, possibly 



because P Gc depends stronger on A^. 

In conclusion, we have considered the Coulomb scat- 
tering between localized orbits as an electronic process 
for the energy dissipation in the breakdown of the quan- 
tum Hall effect in the presence of slowly-fluctuating po- 
tentials. Compared to the previous theories based on 
inter-Landau-level phonon scatterings, we have obtained 
better agreements with the experiment in the value of 
the lower critical electric field, and in the existence of 
the bistability. 
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